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Simple  mixture  models  for  the  prediction  of  shock  equations  of  state  (Hugoniot)  are  a  necessary 
tool  for  characterization  of  multiple  composites.  A  mixture  model  for  determining  the  shock 
equation  of  state  of  composite  materials  is  presented.  The  model  is  completely  flexible  allowing  for 
multiple  (>2)  components.  Additionally,  error  propagation  analysis  for  the  two  component  mixture 
model  has  been  accomplished.  The  model  predicts  the  equation  of  state  to  5% — 15%  of  the 
experimental  data,  which  is  comparable  to  variations  realized  in  meso-scale  modeling  of  similar 
materials.  [http://dx.doi.org/10.1063/L3702873] 


I.  INTRODUCTION 


II.  MIXTURE  MODEL 


As  composites  become  increasingly  used  in  applications 
where  they  are  exposed  to  shock  loading,  a  mixture  model 
able  to  predict  the  shock  equation  of  state  based  on  the  prop¬ 
erties  of  the  constituents  is  necessary. 1-1 1  However,  models 
presented  in  the  literature  primarily  account  for  only  two 
component  composites,1-6  are  restricted  by  geometry,7  or 
require  finite  element  software. 4-6'x  Bennett  and  Horie  pres¬ 
ent  a  detailed  mixture  model  for  reactive  powder  composites, 
including  calculating  the  volumetric  mixture  properties  of 
the  final  compacted  intermetallic  material.14  Current  trends 
in  mesoscale  modeling  have  resulted  in  predictive  capability 
for  the  properties  of  composites, 15-17  which  require  either 
high  quality  microstructures  for  the  materials  or  the  ability  to 
artificially  generate  such  microstractures.  Given  the  differen¬ 
ces  in  materials,  particularly  for  polymer-based  particulate 
composites  consisting  of  hard  particles  in  soft  binders,  the 
preparation  of  samples  for  microstructural  analysis  can  be 
difficult. 

There  is  very  little  discussion  in  the  literature  about  the 
properties  of  the  individual  constituents,  and  how  well  they 
are  known,  affect  the  predicted  properties  of  the  composites. 
Often  the  models  used  are  so  complex  as  to  prohibit  the 
propagation  of  error  from  the  constituent  input  terms. 

This  paper  presents  a  non-reactive  mixture  model  for 
fully  dense  composites,  which  can  be  implemented  in  a 
mathematical  code,  such  as  Mat  lab.1  The  model  is  flexible 
in  the  number  of  components  that  can  be  accommodate  and 
is  non-specific  in  the  types  of  materials,  e.g.,  polymer  matrix 
composites  and  metallic  alloys  can  both  be  accommodated. 
The  error  in  the  constituent  properties  has  been  propagated 
through  the  two-constituent  model  to  determine  the  error 
bars  for  the  calculated  equation  of  state.  Several  two-  and 
three-constituent  composites  are  investigated  with  the  model 
and  compared  with  experimental  data.  The  paper  builds 
upon  previous  work  which  investigated  the  equation  of  state 
of  epoxy-based  particulate  composites.19 


A.  Derivation  of  mixture  model 

In  determining  the  nonreactive  EOS  for  a  multiphase 
mixture,  mixture  theory  has  been  successfully  applied  to  a 
variety  of  multiphase  materials.  A  mixture  shock  Hugoniot 
is  defined  by  determining  the  relationship  of  shock  velocity 
to  particle  velocity,  Us  and  up,  respectively. 

Using  jump  state  and  characteristics  theory,1  the  shock 
velocity  is,  by  definition,  Us  =  cs  +  up  and  across  a  shock 
the  change  in  density,  for  a  two  constituent  mixture,  is  con¬ 
strained  by 


0i V?  +  <p2P°2  1  _ 

0lPl  +  02P2  Us  ’ 


where  (f>k(k=  1 ,2)  is  the  volume  fraction,  pk  (k=  1 ,  2)  is  the 
true  material  density,  and  the  0  superscript  corresponds  to 
the  undisturbed  state.  In  the  limit  of  homogeneous  flow  for  a 
two  constituent  mixture,  the  mixture  sound  speed,  cs,  is 
given  by20 


1 

(Us  -  up)2 


—  (0lPl  +  02P2) 


(2) 


For  a  saturated  mixture. 


01  +  02—1 


(3) 


and  the  true  density  for  each  constituent  is  determined  using 
the  shock  Hugoniot  for  the  pure  constituents, 


P°  _  !  up 
Pk  Us,k 


(4) 


and 


Us,k  —  CQ,k  +  S\Up,k  +  S2Up  t 


(5) 


a,Author  to  whom  correspondence  should  be  addressed.  Electronic  mail: 
Jennifer.jordan@eglin.af.mil. 


Therefore,  a  mixture  shock  velocity  (corresponding  to  a 
specified  particle  velocity)  is  determined  by  combining 
Eqs.  (1)  and  (2)  and  the  resulting  transcendental  relationship 
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(for  a  specific  volume  fraction)  is  solved  using  a  standard 
root  solving  method.  The  mixture  quantities  for  pressure,  P, 
and  density  are  then  defined  by 

P  =  P°  +  (0?p?  +  (j}°2p02)upUs  (6) 


a/1 

1  /?(1  —  a)  +  a  ’ 

/?(1  -  a) 


02  ~ 


/?(1  -  a)  +  a  ’ 


(16a) 

(16b) 


and  Eq.  (1). 

For  a  three  constituent  mixture,  similar  relationships 
hold  and  the  procedure  is  applied  first  treating  two  constitu¬ 
ents  as  a  new  “1  +  2”  constituent  given  that  </>j  >  0;  hence, 
one  defines  a  mixture-averaged  density  and  sound  speed,  p12 
and  c i2-  A  normalized  volume  fraction  is  given  by 


01  +  02 

and 

Pn  =  “Pi  +  (1  -  «)P2  (8) 

the  mixture  relationships  of  “  1  +  2”  phases  are  similar  to 
Eqs.  (2)  and  (4)  rewritten  in  terms  of  the  normalized  volume 
fraction, 


and 


«(1  -  ft) 

/?(1  -  a)  +  a  ’ 


(16c) 


Conservation  of  mass  across  the  jump  state  (shock)  requires 


0Vl  +  02f2  +  03  Pi  1  Up 

0lPl  +  02P2  +  03P3  Us,  123 


(17) 


Hence,  the  wave  speed  for  a  three  phase  mixture,  Us323,  is 
defined  as  a  function  of  particle  velocity,  up,  producing  the 
shock  Hugoniot  for  a  three  phase  system. 

A  similar  procedure  is  applied  to  a  multi-component 
mixture  containing  four  or  more  constituents.  In  this  general 
procedure,  a  reference  component,  subscripted  by  “1”  such 
that  0?>  0,  and  normalized  variables  are  defined. 


P\2  _  j  UP  _  ^12 

Pl2  £4,12  Up  +  Ci2 


(9) 


and 


C12  (£7.5,12  —  Up) 


■=  («Pi  +  (1  -  oc)p2)  ( 


1  —  a\ 


pA 


piA 


(10) 


and 


0' 

01+0/ 

Plj  =  <XjPj+  (1  -  <Xj)Pu 

Py  _  j  _  up  _  cy 
Pij  £7.v. i j  Up  T  c ] j 


(18) 

(19) 

(20) 


Given  the  particle  velocity,  up ,  true  densities  for  1  and  2  are 
defined  by  their  respective  EOS  in  their  pure  state  and  a  is 
solved  combining  Eqs.  (9)  and  (10)  and  solving  the  resulting 
equation  with  a  root  solver. 

Similar  relationships  are  used  for  a  new  “1  +  3”  constit¬ 
uent  defining  p13  and  c13  by  the  following  equations: 


P  = 


0i  +  03  ' 

Pl3  =  PPl  +  (1  -  /?)P3, 
U„  C 13 


P 13  _  ^  _  "P 
Pl3 


(ID 

(12) 

(13) 


1 

X 

-13 


£7S,  13  Up  +  C 13  ’ 

(JJ— - ,2  =  f/^Pl  +  (1  -  P)P3)  (A  +  1— I) 

(C7,,13  —  Mp)“  V  J\P\C\  PiCiJ 


(14) 


1 

(£/S,ly  -  Up)~ 


(21) 


The  solution  for  all  otj  (where  y.j  <  I )  can  be  determined  upon 
specifying  equations  of  state  for  each  individual  constituent 
and  solving  the  sets  of  transcendental  Equations  (18) — (21). 
Since  saturation  implies  that  ,  <f)k  =  1 ,  then 

0*  =  ^^,  for  A- =  2,  TV  (22) 

1  3-k 


hence. 


1 


(23a) 


and  (1  is  solved  combining  Eqs.  (13)  and  (14),  and  the  result¬ 
ing  equation  is  solved  using  a  root  solver. 

Finally,  for  a  three  constituent  system,  the  saturation 
constraint  requires 


0i  +  02  +  03  —  1  (15) 

and  with  Eqs.  (7)  and  (11),  the  shocked  volume  fractions  are 
defined  for  a  specified  Mp 


and 


(1  °v  )  (  1  + 


£2 1  -  **, 


(23b) 


and  the  overall  density  change  is  then  related  to  the  mixture 
wave  speed  to  define  a  mixture  equation  of  state  (i.e.,  shock 
Hugoniot)  according  to 
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k=l 


N 

E  <PkPk 

k=  1 


1  -- 


£4. 


(24) 


B.  Uncertainty  analysis  in  mixture  relationships 

Along  the  Hugoniot,  the  mixture  pressure,  P,  is  related 
to  the  average  particle  velocity,  up,  and  the  mixture  wave 
speed,  Us,  by 


fttupU, 


(25) 


(U]-Usup)2(  -(1  ~$)p°2h  \ 

Up  ~  2US  \(p®)~Us\  (17.51  _  Up)  P\Usl{Us2  —  Up)  J 

(31b) 

W-usup)2(  (1-^)0,  \ 

up  ~  2^4  \Pi Usi(Usi  —  Up)  (p2)~Us2(Us2  —  up) ) 

(31c) 


( U)  ~  UsUpf  (hPo\  (  2 UA-Up  \ 

up~2Us  V  p\  )\{Ul-Uslup)2)' 


(3  Id) 


where  p™lx  is  the  undisturbed  density.  At  a  fixed  particle  ve¬ 
locity,  it  is  readily  seen  that 

A P  A  Us 

—  = - -  (26) 

P  Us 

or  the  uncertainty  associated  with  estimating  pressure  is 
directly  related  to  the  uncertainty  in  the  mixture  wave  speed. 
Hence,  the  variations  of  mixture  wave  speed  are  assessed 
with  regard  to  uncertainties  of  the  individual  constituents  of 
the  mixture. 

For  a  two  component  mixture,  consider 


Us=f(p°l,p02,(l>0l,(l>1,UshUS2) 


(27) 


[U]  -  Usupf  /(1-(Mpq\  (  2U*i~up  \ 

“P-  2US  \  p°2  )  {(U^2  -  Us2up)2  J  ’ 

(31e) 


a  5  = 


(U2S  -  Usupf 
Up  2  Us 


(pI-pI) 


(  01  ,  1-01 
\AU,i(Usi  -  Up)  p°2us2(us2  -  Up) 


(3  If) 


a6 


(U's  —  UsUpj2  (  1 

Up  -  2 Us  Po  Xp^UsiiUsi  -  Up) 


P°2Us2(Us2-up)J' 

(3  lg) 


at  conditions  of  fixed  up.  Following  Ref.  21,  the  appropriate 
relationship  for  the  uncertainty  of  mixture  wave  speed,  dUs, 
in  terms  of  the  uncertainties  of  variables  in  Eq.  (27)  is 


dU,= 


+ 


(S». 


dus 

dUs 


-dUsl 


dU, 


dU, 


-dUs  2 


s2 


(28) 


The  function  for  Us  is  defined  by  the  relationship  for  the 
mixture  wave  speed,  i.e.,  Eq.  (2),  which,  knowing  that 


All  of  the  terms  in  Eq.  (31a)  are  known,  except  for  the  error 
in  the  volume  fraction  at  pressure,  d<fil,  since  the  volume 
fraction  at  pressure  is  a  derived  quantity.  However,  this  error 
can  be  determined  by  differentiating  Eq.  (1), 


d(f>l  =  {(bidUs)2  +  (b2d(j)[ l)2  +  {b3dp°\)2  +  (b4dp2)2 

2  9't1/2 

+  (b5dUsi)  +  {b6dUs2)  j  .  (32a) 


where 


b\ 


(32b) 


P  _  Us 
Po  U s  up 


(29) 


bi_  (P\-P\)U* 

(£4  Up )  fp)bp  1 2 


(32c) 


can  be  rearranged  as 

1  =  ^  ] 
US(US  —  up)  \ p°Usi(Usi  -  up)  p\U s2{U s2  Up)  J  " 

(30) 

Partial  differentiation  of  Eq.  (30)  in  terms  of  each  of  the  vari¬ 
ables  in  Eq.  (28)  yields 

dUs  +  ( a2dp\ )2  +  (a3dUsi)2 

+  (a4dUs2)2  +  (a2d<j)0l)2  +  (a6<70i)2)  ^  ,  (31a) 

where 


b3 


b  4 


PqUs  _  P°2Us 2  \  /  Us  1  \ 

44  Up  Us  2  Up  J  \  ( U s  ]  Up)  J 

-(p?)2 


PoUs  _  p\UsX  4  /  Us2  \ 

(74  Up)  (74i  Up) J  \(Us2  Up) J 

-{bpl2Po)2 


(32d) 


(32e) 


bs 


(  PqUs  P°2Us2  \ 

|  pf\llp  \  \{U s  ~  Up)  (U s2  Up)) 

V(£4t  -  Up)2)  -(bpnPof 


(32f) 
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TABLE  I.  Constituent  and  composite  parameters  for  uranium  with  1 3.4  wt.  % 
rhodium. 


Material 

Po 

(g/cm3) 

C0 

(km/s) 

S 

0° 

dpo 

d(j>0 

dUs 

Refs. 

U 

18.93 

2.51 

1.51 

0.735 

0.05 

0.001 

0.061 

22 

Rh 

12.429 

4.79 

1.36 

0.265 

0.05 

0.058 

22 

U-13.4  wt.  %  Rh 

17.20 

2.74 

1.49 

22 

b6  = 


PlllP 


PqUs 


P°iUs 


XU  a  -  uPY 


( Us  Up )  ( Us  i  Up) 

(bpuPo)2 


(32g) 


and 


dpi2  — 


P  0 


P°lUsl  P°2Us  2 

(C/,1  Rp)  ( Uv2  Up ) 


(33) 


Equation  (30)  can  be  substituted  into  Eq.  (31a)  and  rear¬ 
ranged  to  solve  for  dUs, 


dUs  = 


1  -  alh  i 

,2.2 


Ot  +  «6^3)(^i)2  +  («2  +  a&Didp'i) 


,0\2 


+  (a3  +  a6^55)  Wl )  +  (fl4  +  a6^6)  s2)“ 


+  (ag  +  a\bY)(d<f>\) 


/  0\2 


1/2 


(34) 


The  error  propagation  for  systems  containing  three  or  more 
components  is  increasingly  complex  and  is  not  considered  in 
this  paper. 


III.  COMPARISON  WITH  EXPERIMENTAL  DATA 
A.  Two-constituent  composites 
1.  Uranium-rhodium  composite 

Uranium-rhodium  (U-13.4wt.  %  Rh)  was  chosen  as  a 
model  metal-metal  composite.22  The  properties  for  each  con¬ 
stituent  used  in  the  model  are  given  in  Table  I.  The  bulk 
sound  speed,  C0,  and  S  were  determined  from  linear  least 
squares  fit  to  the  available  data22  as  shown  in  Figs.  1(a)  and 
1(b)  for  uranium  and  rhodium,  respectively.  The  deviation 
from  these  linear  fits,  dUs,  were  also  determined  and  are 
shown  in  Fig.  1  and  Table  I. 

The  predicted  Hugoniot  relationship  for  U- 13.4  wt.  %  Rh 
was  determined  from  the  mixture  model  for  two  constituents 
presented  above.  It  can  be  seen  from  Fig.  2(a)  that  the  model 
slightly  overpredicts  the  experimental  data,  with  an  average 
deviation,  dUJUs  of  0.05,  shown  in  Fig.  2(b).  The  linear  fits 
for  uranium  and  rhodium  are  shown  for  comparison. 

The  propagated  error  associated  with  the  input  parame¬ 
ters  to  the  mixture  model  was  determined  from  Eq.  (34)  and 
is  shown  as  the  dashed  lines  bracketing  the  solid  line  in  Fig. 
2(a).  The  error  in  density  was  determined  from  the  differen¬ 
ces  in  the  reported  density  values  for  uranium,  and  the  same 
error  was  used  for  the  density  in  rhodium.  The  density  of  the 
composite  determined  from  the  mixture  model  calculation 
agrees  with  that  reported  for  the  composite.22  In  order  to 
determine  the  error  in  the  input  volume  fraction,  it  was 
assumed  that  there  was  +/—  1  unit  error  in  the  last  reported 
significant  figure.  The  error  in  the  shock  velocity  for  both 
constituents  is  the  dUs  values  shown  in  Figs.  1(c)  and  1(d). 
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FIG.  1 .  Particle  velocity  vs  shock  veloc¬ 
ity  for  (a)  uranium  and  (b)  rhodium, 
where  the  circles  are  the  experimental 
data  from  Ref.  20  and  the  solid  line  is  a 
linear  fit  to  these  experimental  data,  and 
dUs/Us  for  (c)  uranium  and  (d) 
rhodium. 
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FIG.  2.  (a)  Particle  velocity  vs  shock 
velocity  and  (b)  dUs/Us  for  uranium, 
13.4  wt.%  rhodium  experimental  data 
(Ref.  20)  compared  with  the  mixture 
model  calculation.  The  dashed  lines  on 
either  side  of  the  solid  line  in  (a)  indi¬ 
cated  the  error  bars  for  the  model 
calculation. 


The  experimental  data seem  to  lie  on  or  around  the  lower 
error  bar  for  this  composite. 

2.  Epoxy-aluminum  and  epoxy-alumina  composites 

Two  epoxy-based  two-phase  particulate  composites 
were  considered:  Al203-57  vol  %  epoxy23-26  and  Al-55  vol 
%  epoxy.27  The  former  material  has  been  extensively  studied 
in  the  literature  and  the  experimental  data  points  used  are 
combined  from  multiple  references.23-26  The  latter  material 
had  a  limited  number  of  experimental  Hugoniot  points  avail¬ 
able,  but  a  comparison  between  a  polymer-based  composite 
containing  a  metal  powder  and  one  containing  an  oxide  pow¬ 
der  should  be  useful. 

The  properties  of  the  constituent  materials  are  presented 
in  Table  II  (Refs.  22-32)  and  Fig.  3.  The  epoxy  used  in  this 
study  appears  to  undergo  a  phase  change  around 
up  =  2.84  km/s,  Fig.  3(a).  In  order  accommodate  this  phase 
change  as  well  as  a  downturn  in  Us  at  low  values  of  up,  a 
piecewise  fit  to  the  Us  -  up  data  for  epoxy  were  used  and  is 
listed  in  Table  III.  In  order  to  propagate  the  error  from  the 
input  parameters  through  the  model  to  determine  the  error 
bars  in  the  output,  the  variations  in  the  input  parameters 
must  be  known.  The  variation  in  shock  velocity,  dUs,  was 
determined  from  the  difference  between  the  experimental 
data  points  and  the  linear  least  squares  fit  to  the  Us  -  up  data. 
This  variation  is  presented  graphically  in  Figs.  3(d)— 3(f)  and 


tabulated  in  Table  II.  Alumina  deviates  significantly  from  a 
linear  Us—up  fit  in  the  low  pressure  regime.  Here  the  average 
dUs,  used  in  propagating  the  error  for  the  mixture,  may  not 
be  representative.  The  variation  in  density  for  each  of  the 
materials  was  determined  from  either  the  variation  in 
reported  densities  in  the  literature  or  assuming  a  variation  of 
+/— 1  unit  in  the  last  reported  significant  figure. 

The  model  prediction  along  with  the  experimental  data 
for  AI2O3-57  vol  %  epoxy  is  presented  in  Fig.  4(a)  and  for 
Al-55  vol  %  epoxy  in  Fig.  4(c),  along  with  the  linear  fits  for 
the  constituents.  The  model  prediction  underpredicts  the 
data  for  Al203-57  vol  %  epoxy  by  approximately  14% 
(Fig.  4(b)).  The  prediction  for  Al-55  vol  %  epoxy  is  better  at 
approximately  5%;  however,  there  exists  limited  data  points 
available  to  compare  to  the  model  predictions. 

B.  Multi-component  composites 
1.  304  stainless  steel 

304  stainless  steel  was  chosen  as  a  representative  three 
constituent  composite,  with  the  primary  constituents,  74% 
Fe,  18%  Ni,  and  8%  Cr,33  considered  in  this  study.  The 
phase  change  in  iron  was  neglected  in  the  model,  although  it 
could  be  incorporated  in  the  same  method  as  the  epoxy 
above.  The  properties  for  each  constituent  used  in  the  model 
are  given  in  Table  IV.  The  bulk  sound  speed,  C0,  and  5  were 


TABLE  II.  Constituent  and  composite  parameters  for  epoxy-aluminum  and  epoxy-alumina  materials. 


Material 

p0  (g/cm3) 

C0  (km/s) 

S 

$ epoxy 

dpo 

d<t>° 

dUs 

Ref. 

Epoxy 

1.134—1.194“ 

0.03 

0.09 

28-32 

Aluminum 

2.712 

5.33 

1.38 

0.001 

0.06 

22 

Alumina 

3.976 

8.14 

1.28 

0.007 

0.24 

22 

Epoxy-Al203  (300/100/20) 

2.37 

2.93 

1.94 

0.58 

0.01 

23 

Epoxy-Al203  (A12) 

2.377 

2.88 

1.99 

0.57 

0.01 

24 

Epoxy-Al203  (A13) 

2.389 

0.57 

0.01 

24 

Epoxy-Al203  (A14) 

2.391 

0.57 

0.01 

24 

Epoxy- A1203  (fully  loaded) 

2.28 

2.93 

1.63 

0.57 

0.01 

25 

Epoxy- A1203  (FI) 

2.376 

2.76 

2.26 

0.01 

26 

Epoxy- AI2O3  (F2) 

2.429 

2.02 

2.26 

0.01 

26 

Epoxy- A1  (Epoxy-45X8 1 ) 

1.875 

0.55 

0.01 

27 

aDependent  on  resin  and  curing  agent. 
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TABLE  III.  Input  parameters  for  mixture  and  mesoscale  models. 


Parameter 

Epoxy 

Phase  I 

Epoxy 

Phase  II 

Epoxy 

Phase  III 

Epoxy 

Phase  IV 

Up  range 

0-0.4 

0.4-2.81 

2.81-3.1 

3. 1-5.2 

C0  (km/s) 

2.263 

2.67 

7.17 

2.88 

S 

3.4704 

1.55 

0.01 

1.35 

S2 

-0.3947 

determined  from  linear  least  squares  fits  for  Ni  and  Cr  (Ref. 
22)  and  used  as  reported  in  Ref.  34  for  Fe. 

The  predicted  Hugoniot  relationship  for  304  stainless 
steel  was  determined  from  the  mixture  model  for  three  con¬ 
stituents  presented  above.  It  can  be  seen  from  Fig.  5  that  the 
model  slightly  underpredicts  the  experimental  data,  with  an 
average  deviation,  dUs,/Us  of  0.05,  shown  in  Fig.  5(b).  The 
linear  fits  for  iron,  nickel,  and  chromium  are  shown  for  com¬ 
parison.  The  model  may  be  improved  by  the  addition  of  the 
minor  components  in  304  stainless  steel  and  inclusion  of  the 
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FIG.  3.  Particle  velocity  vs  shock  veloc¬ 
ity  for  (a)  epoxy,  (b)  alumina,  and  (c) 
aluminum,  where  the  circles  are  the  ex¬ 
perimental  data  from  Refs.  20,  26-30, 
and  the  solid  line  is  a  linear  least  squares 
fit  to  this  experimental  data,  and  dUs/Us 
for  (d)  epoxy,  (e)  alumina,  and  (f) 
aluminum. 
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FIG.  4.  (a)  Particle  velocity  vs  shock  ve¬ 
locity  and  (b)  dUs/Us  for  Al203-57  vol 
%  epoxy  and  (c)  particle  velocity  vs 
shock  velocity  and  (d)  dUs/Us  for  Al-55 
vol  %  epoxy  experimental  data  com¬ 
pared  with  mixture  model  calculation. 
The  dashed  lines  on  either  side  of  the 
solid  line  in  (a)  indicated  the  error  bars 
for  the  model  calculation. 


phase  change  in  iron;  however,  the  agreement  for  this  three 
component  system  is  as  good  as  that  for  the  uranium- 
rhodium  composite. 

2.  Epoxy-Al-Fe203  and  epoxy-AI-Mn02 

Two  multi-component  epoxy-based  composites  have 
sufficient  experimental  data  available  in  the  literature,  ep- 
oxy-Al-Fe2C>3  (Ref.  35)  and  epoxy-Al-Mn0276  The  input 
data  for  epoxy  and  aluminum  have  been  discussed  in  the  pre¬ 
vious  section.  The  equation  of  state  for  Mn02  was  deter¬ 
mined  from  that  for  pyrolusite,  as  presented  in  Ref.  22.  All 
of  the  input  data  for  the  multi-phase  composites  is  given  in 
Table  V.  The  equations  of  state  for  Fc202  and  Mn02  are  pre¬ 
sented  graphically  in  Fig.  6  with  their  associated  variation  in 
shock  wave  speed.  For  both  materials,  epoxy  was  chosen  as 
material  “1,”  or  the  reference  material,  for  Eqs.  (7 ) — ( 15). 
However,  the  model  is  insensitive  to  the  choice  of  the  refer¬ 
ence  material,  and  simulations  of  the  epoxy-Al-Fe203  system 


TABLE  IV.  Constituent  and  composite  parameters  for  304  stainless  steel 
and  its  constituents. 


Material 

p0  (g/cm3) 

C0  (km/s) 

S 

Refs. 

Fe 

7.85 

3.57 

1.92 

0.74 

34 

Ni 

8.875 

4.59 

1.44 

0.18 

22 

Cr 

7.119 

5.2 

1.43 

0.08 

22 

304  SS 

7.89 

4.56 

1.50 

22 

were  conducted  with  each  constituent  as  the  reference  mate¬ 
rial  with  no  change  in  the  predicted  mixture  behavior. 

The  experimental  data  and  predicted  model  for  epoxy- 
Al-Fe203  is  shown  in  Fig.  7(a)  along  with  the  linear  fits  to 
the  equations  of  state  for  the  constituents.  It  can  be  seen  that 
there  is  reasonably  good  agreement  between  the  model  and 
the  experimental  data.  The  ~7%  error.  Fig.  7(b),  may  be  pri¬ 
marily  due  to  the  scatter  in  the  experimental  data  rather  than 
systematic  differences  between  the  experiments  and  the  sim¬ 
ulation.  The  model  and  experimental  data  are  very  close  to 
the  equation  of  state  for  epoxy,  which  is  most  likely  due  to 
the  high  volume  fraction  of  epoxy  in  this  system. 

The  predicted  model  and  the  experimental  data  for  ep- 
oxy-Al-Mn02  are  presented  in  Fig.  7(c).  For  this  composite, 
the  density  from  the  mixture  model  did  not  have  good  agree¬ 
ment  with  the  actual  density  measured  for  the  composite. 
The  Mn02  powders  displayed  porosity  in  Ref.  10,  so  the  ini¬ 
tial  density  of  the  Mn02  constituent  was  modified  by  match¬ 
ing  the  composite  mixture  density  to  the  experimentally 
measured  composite  density.  The  adjusted  density  is  shown 
in  parentheses  in  Table  IV.  With  the  adjusted  density  the 
agreement  between  the  experimental  data  and  the  model  is 
approximately  13%,  which  is  in  line  with  the  agreement  with 
two-phase  epoxy  based  composites. 

The  experimental  data  in  this  system  lies  above  the  ep¬ 
oxy  equation  of  state,  rather  than  around  it  as  in  the  epoxy- 
Al-Fe203  system.  This  indicates  that  the  decreased  volume 
fraction  of  epoxy  in  this  system  is  reducing  the  effect  of  the 
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FIG.  5.  (a)  Particle  velocity  vs  shock  velocity 
and  (b)  dUs/Us  for  304  stainless  steel,  experi¬ 
mental  data  (Ref.  20)  compared  with  the  mixture 
model  calculation. 


FIG.  6.  Particle  velocity  vs  shock  velocity  for 
(a)  Fe203  and  (b)  Mn02,  where  the  circles  are 
the  experimental  data  from  Ref.  20  and  the 
solid  line  is  a  linear  fit  to  these  experimental 
data,  and  dUs/Us  for  (c)  Fe203  and  (d) 
Mn02. 


TABLE  V.  Constituent  and  composite  parameters  for  epoxy-Al-Fe203  and  epoxy-Al-Mn02  materials. 


Material 

Pa  (g/cm3) 

C0  (km/s) 

S 

/  0 

t epoxy 

dpo 

d<j>° 

dUs 

Refs. 

Epoxy 

1.134—1.194“ 

0.03 

0.09 

28-32 

Aluminum 

2.712 

5.33 

1.38 

0.001 

0.06 

22 

Iron  oxide  (Fe203) 

5.00 

6.24 

1.39 

0.01 

0.21 

22 

Manganese  dioxide  (Mn02) 

5.03  (4.55) 

3.261 

1.24 

0.01 

22 

Epoxy-Al-Fe203 

2.366 

0.60 

0.01 

35 

Epoxy-Al-Mn02 

2.60 

3.61  b 

1.47  b 

0.45 

0.02 

0.01 

36 

aDependent  on  resin  and  curing  agent. 
bLinear  fit  0.7  <up<  1.7  km/s. 
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FIG.  7.  (a)  Particle  velocity  vs  shock  ve¬ 
locity  and  (b)  dUs/Us  for  epoxy-Al- 
Fe203  and  (c)  particle  velocity  vs  shock 
velocity  and  (d)  dUs/Us  for  epoxy-Al- 
M11O2  experimental  (Refs.  30  and  31) 
data  compared  with  the  mixture  model 
calculation.  The  dashed  lines  on  either 
side  of  the  solid  line  in  (a)  indicated  the 
error  bars  for  the  model  calculation. 


epoxy  binder  on  the  composite  behavior.  The  model  does  not 
reflect  this  same  change  in  the  influence  of  the  epoxy  binder. 
Increasing  the  aluminum  volume  fraction  in  the  model  to 
0.40  (versus  0.24  actually  in  the  composite),  and  decreasing 
the  epoxy  volume  fraction  by  a  comparable  amount,  results 
in  a  model  that  more  accurately  reflects  the  behavior  of  the 
composite.  This  may  indicate  that  the  volume  fraction  of  the 
components  varied  in  the  samples  or  that  the  equation  of 
state  approximated  for  MnCT  is  not  accurate. 

IV.  CONCLUSIONS 

A  general  mixture  model  for  determining  the  shock 
equation  of  state  of  composite  materials  is  presented.  The 
model  is  completely  flexible  allowing  for  multiple  (>2)  con¬ 
stituents  in  the  composite.  Additionally,  error  propagation 
analysis  for  the  two  component  mixture  model  has  been 
accomplished.  For  metallic  composites,  i.e.,  uranium- 
rhodium  and  304  stainless  steel,  the  average  error  between 
experiment  and  model  is  ~5%.  The  prediction  error  for 
epoxy-based  composites  increases  to  ^10%-15%  based  on 
the  uncertainty  in  the  epoxy  equation  of  state,  which  domi¬ 
nates  the  material  behavior  due  to  its  relatively  large  volume 
fraction.  This  mixture  model  allows  for  prediction  of  the 
equation  of  state  of  complex  composite  materials  within 
5%-15%  of  the  experimental  data,  which  is  consistent  with 
variations  that  are  realized  in  meso-scale  modeling  of  similar 
materials.  The  model  can  be  implemented  in  MatLab  and 


run  in  minutes  on  a  desktop  computer,  which  may  provide 
an  advantage  over  the  longer  running  meso-scale  simula¬ 
tions,  which  also  require  significant  computational  set-up  to 
be  able  to  run.  This  mixture  model  tool  may  be  used  to 
predict  experiments  involving  complex  materials  that  have 
not  been  previously  tested,  allowing  for  confidence  in  experi¬ 
mental  timing  and  other  parameters. 
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